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ABSTRACT 

We present observations of e Eridani from the Submillimeter Array (SMA) 
at 1.3 millimeters and from the Australia Telescope Compact Array (ATCA) 
at 7 millimeters that reach an angular resolution of ~ 4" (13 AU). These hrst 
millimeter interferometer observations of e Eridani, which hosts the closest de¬ 
bris disk to the Sun, reveal two distinct emission components: (1) the well-known 
outer dust belt, which, although patchy, is clearly resolved in the radial direction, 
and (2) an unresolved source coincident with the position of the star. We use 
direct model-htting of the millimeter visibilities to constrain the basic properties 
of these two components. A simple Gaussian shape for the outer belt £t to the 
SMA data results in a radial location of 64.4^(3 q AU and FWHM of 20.2^g 2 AU 


2 


(fractional width AR/R = 0.3). Similar results are obtained taking a power law 
radial emission profile for the belt, though the power law index cannot be usefully 
constrained. Within the noise obtained (0.2 mJy beam“^), these data are consis¬ 
tent with an axisymmetric belt model and show no signihcant azimuthal structure 
that might be introduced by unseen planets in the system. These data also limit 
any stellocentric offset of the belt to < 9 AU, which disfavors the presence of 
giant planets on highly eccentric (> 0.1) and wide (lO’s of AU) orbits. The flux 
density of the unresolved central component exceeds predictions for the stellar 
photosphere at these long wavelengths, by a marginally signihcant amount at 
1.3 millimeters but by a factor of a few at 7 millimeters (with brightness temper¬ 
ature 13000 ± 1600 K for a source size of the optical stellar radius). We attribute 
this excess emission to ionized plasma from a stellar corona or chromosphere. 

Subject headings: circumstellar matter — stars: individual (e Eridani) — sub- 
millimeter: planetary systems 


Introduction 


Debris disks, composed of planetesimals remaining after planet formation and circum- 
ste llar disk dispersion, represent the end-stage of protoplanetary di sk evolution (see reviews 


by iRackman fc Parescd Il993l: IWvattI l2008l: Matthews et al.l 120141) . While these remnant 


planetesimals cannot be observed directly, they are ground down through ongoing collisions 
into smaller and smaller dust grains that scatter starlight and produce detectable thermal 
emission. Observations of this dusty debris at millimeter wavelengths are especially critical 
to our understanding of the most readily accessible systems. The large grains that dominate 
emission at these long wavelengths do not travel far from their origin and therefore reliably 
trace the underlying planetesimals di stribution, un like the small grains that are rapidly re¬ 


moved by stellar radiation and winds fjWvattll2006[) . Since planets, if present, will inevitably 


perturb the dust-producing planetesimals, the millimeter emission morphology encodes in¬ 
formation on the architecture and dynamical evolution of these systems. For example, the 
outward migration of a p lanet can conhne planetesimals in a belt between its resonances 


orbit ( 

Kuchner & Holman 

2003: 

Wvatt 

2003; 

Deller & Maddison 

2005) 

sculpt out sharp edges in a belt 

f Onille 

n boOl 

il: Chiane et ah 

200 

9), or 

onto eccentric or inclined orbits i 

Wvatt et ah 

1999 

)• 


At a distance of only 3.22 pc (Ivan Leeuwenll2007l) . the 400—800 Myr-old flMamaiek fc Hillenbrand 


20081) main-sequence K2 star e Eridani hosts the closest debris disk to the Sun, originally 
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identified through the detection of far-infrared emission by IRAS flAnmannlh985h . Pioneering 
observations with JCMT/SCUBA resolved a nearly face-on belt of emission at 850 /rm, peak¬ 


ing a t 60 AU (18") radius, with several brightness enhancements or clumps flGreaves et ah 


19981) . Analysis bv iGreaves et all (120051) of JGMT images spanning 5 years offered tentative 


evidence that some of these clumps are stationary, and so likely background galaxies, while 
others appear to be co-moving with the star, and so are likely associated with the disk. 
Subsequent single-dish observations from 250 to 1200 pm have conhrmed the basic belt mor¬ 


phology, but not all of the low signihcance asymmet ries (jSchutz et al.ll2004J: iBackman et ah 


20091: iGreaves et al.ll2014J: iLestrade &: ThilliezI 120151) . A warmer dust component, reaching 


to several AU from the star , can be explained in modeling the spectr al energy distribution 


by an additional dust belt (IBackman et al 


2009 


Greaves et ahl 120141) . or by inward trans¬ 


port from the outer belt (iReidemeister et al.N201l[) . In addition, precision radial velocity 
obser vations suggest the presence of a Jupiter-mass planet with semi-major axis of 3.4 AU 
(1") (IHatzes et al.ll2000l). althou gh the reality of this planet signal remains controversial 
(lAnglada-Escude fc Butleiil2012l) . Because e Eridani is so nearby, it is a key template for 
understanding debris disk phenomena around Sun-like stars, and detailed study of its debris 
disk provides essential context for the interpretation of more distant, less accessible systems. 


We present observations of e Eridani at 1.3 mm and at 7 mm, using the Submillimeter 
Array (SMA) and the Australia Telescope Gompact Array (ATGA), respectively, using the 
most compact, lowest angular resolution (4" — 10") conhgurations of these telescopes. While 
even higher resolution may be desirable, these interferometric observations are conservatively 
tuned to provide a hrst look at structures below the resolution of previous single dish obser¬ 
vations. For these arrays at these wavelengths, the primary beam held of view encompasses 
the entire emission region from the outer debris belt, enabling efhcient observations of the 
full disk with a single pointing. Section |2] describes these observations of the e Eridani sys¬ 
tem. Section [3] describes the modeling procedure and the results. Section 0] discusses the 
implications of the model hts for the outer dust belt properties, azimuthal asymmetries, and 
the nature of an inner component of excess millimeter emission. 


2. Observations 


2.1. Submillimeter Array 


We observed e Eridani in July, August and November 2014 with the SMA (IHo et al. 


20041 ) on Manna Kea, Hawaii at a wavelength of 1.3 mm in the subcompact conhguration. 
Table [U summarizes the essentials of these observations, including the dates, baseline lengths, 
and atmospheric opacity. Six tracks were obtained, all with 7 operational antennas in the 
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array. The weather conditions were very good for observations at this wavelength (225 GHz 
opacities from 0.07 to 0.12). The total bandwidth available was 8 GHz consisting of two 
sidebands of 4 GHz width spanning ±4 to 8 GHz from the local oscillator (LO) freqnency 
of 225.5 GHz (217.5 — 221.5 GHz and 229.5 — 233.5 GHz). The phase center was located at 
a = 03^32™54!9024, b = —09°27'29'.'4486 (J2000), correspo nding to the positio n of the star 
corrected for its proper motion of (-975.17,19.49) mas yr“^ (Ivan LeeuwenI 120071) as of Jnly 1, 
2014 . At the LO freqnency, the held of view is ~ 52", set by the primary beam size of the 
6 -m diameter array antennas. 


The data from each track were calibrated independently using the IDL-based MIR soft¬ 
ware package. Time-dependent complex gains were determined from observations of two 
nearby quasars, J0339-017 (7?9 away) and J0423-013 (14?9 away), interleaved with observa¬ 
tions of e Eridani in a 16 minute cycle. The passband shape was calibrated using available 
bright sources, mainly 3G84 or 3G454.3. Observations of Uranus or Gallisto during each track 
were used to derive the absolute hux scale with an estimated accuracy of ~ 10%. Imaging 
and deconvolution were performed with the clean task in the GASA software package. A 
variety of visibility weighting schemes were used to explore compromises in imaging between 
higher angular resolution and better surface brightness sensitivity. With natural weighting, 
the beam size is 670 x 575 (19 x 18 AU) and the rms noise level is 0.17 mJy beam“^. The 
longest baselines in the dataset probe size scales of ~ 4" (13 AU). 


2.2. Australia Telescope Compact Array 


We observed e Eridani in late June and early August 2014 with the ATGA, located near 
Narrabri, NSW, at a wavelength of 7 mm using the compact H75 and HI68 conhgurations 
of the array. Table [2] summarizes the essentials of these observations. Four tracks were 
obtained in each of the two antenna conhgurations, all with 6 operational antennas. The 
winter weather provided good atmospheric phase stability for this ATGA high frequency 
band (rms path typically < 150 microns from the seeing monitor), especially for the short 
baselines of interest, except for the last track in the H168 conhguration (2014 July 28), which 
proved to be unusable due to high winds and relatively poor seeing. Data from the stationary 
sixth antenna of the array, located ~ 6 km from the others, was discarded, given the large 
gap from the rest of the antennas and less stable phase on the much longer baselines. The 
bandwidth provided by the Gompact Array Broadband Backend was 8 GHz, with 2 GH z 
wide bands centered at 43 GHz and at 45 GHz, in two polarizations flWilson et al.ll201ll) . 
The phase center was identical to the contemporaneous SMA observations. The held of view 
is ~ 70", set by the primary beam size of the 22-m diameter array antennas. 
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The data from the seven usable tracks were calibrated independently using the miriad 
software package. Time-dependent complex gains were determined using the nearby quasar 
0336-019 (7?8 away), interleaved with observations of e Eridani in a 12 minute cycle. The 
passband shape was calibrated using the available bright sources 1921-293 or 0537-441. Ob¬ 
servations of 1934-638 and Uranus during each track were used to derive the absolute flux 
scale. Comparison of the derived fluxes for 0336-019 for each of the seven nights shows a 
maximum difference of 6%, and we conservatively estimate the flux scale accuracy is better 
than 10%. Imaging and deconvolution were performed with the standard routines invert, 
clean, and restor in the miriad software package. 


3. Results and Analysis 
3.1. Continnum Emission 


Figure [T] shows an SMA 1.3 mm image and ATCA 7 mm image of e Eridani, together with 
the Herschel/SPIRE 250 /im image extracted from the Herschel Science Archive for reference 
flGreaves et al.N20141) . For the 1.3 mm image, the synthesized beam size, obtained with 
natural weighting and a modest taper to improve surface brightness sensitivity is 972 x 877 
(30 X 28 AU), position angle 68°. The rms noise is 0.20 mJy beam“^. This image reveals 
emission from a compact central source at the stellar position (~ 7a) together with patchy 
emission from the nearly face-on belt of cold dust located ~ 18" from the star. For the 
7 mm image, the synthesized beam size obtained with natural weighting is 972 x 770 (30 x 23 
AU), position angle 83°. The rms noise is 7 /iJy beam“^. A central peak is clearly detected 
(~ lOcr). Unlike the 1.3 mm image, the 7 mm image shows little sign of emission from 
the outer dust belt. In both images, the position of the central peak is consistent with the 
predicted stellar position, within the uncertainty dictated by the synthesized b eam size, 9, 


and the signal-to-noise ratio, SNR of ~ QR9/SNR 0'.'6 (see lReid et al.lll988l) . 


3.2. Emission Modeling Procedure 


To characterize the 1.3 m m emission from e Eridani, we used the procedure described by 


MacGregor et al.l (I2013l.l2015h that employs a Markov Chain Monte Carlo (MCMC) method 


to ht simple parametric models to the observed visibilities. We £t the visibility data directly 
both to avoid the non-linear effects of deconvolution and to take advantage of the full range 
of spatial frequencies in the observations that are not necessarily represented well in the 
images. We assume that the emission arises from a geometrically thin, axisymmetric belt. 
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Fig. 1. — fleft) Herschel/SP IRE 250 image of e Eridani, from the Herschel Science Archive 
(see also lGreaves et al.ll2ni4l ). The ellipse in the lower left indicates the ~ 19" beam size, (center) 
SMA image of the 1.3 millimeter continuum emission from e Eridani. The contour levels are in 
steps of 2 X 0.2 mJy, the rms noise level. The ellipse in the lower left corner indicates the T.'2 x 8"7 
(FWHM) synthesized beam size. The dashed blue circle indicates the ~ 52" SMA primary beam 
(FWHM) at the LO frequency (225.5 GHz), (right) ATCA image of the 7 millimeter continuum 
emission from e Eridani. The contour levels are in steps of 2 x 7 /rJy, the rms noise level. The 
ellipse in the lower left corner indicates the 9'.'2 x 7"0 (FWHM) synthesized beam size. The dashed 
blue circle indicates the ~ 70" ATCA primary beam (FWHM) at 44 GHz. In all three panels, the 
star symbol marks the position of the stellar photosphere, corrected for proper motion. For the 
center and right panels, the stellar position is: a = 03^32™54!9024, 6 = —09°27'29'.'4486 (J2000). 


where we consider two different parametric shapes for the snrface brightness profile, Ii,{r): 
(1) an annulus with < r < i?out and power law slope, Rir) oc and (2) a Gaussian, 
R{r) oc exp[—((r — i?cen)/\/2o')^], where Rcen is the position of the belt, a is the width, and 
the FWHM = AR = 2-\/21n(2) x a. The limited signal-to-noise of the dataset precludes 
exploring more complicated, but physically plausible, surface brightness shapes, such as 
multiple rings or a broken power law. The belt emission normalization is defined by a total 
flux density, Fbeit = / Rdfl, and the belt center is given by offsets relative to the pointing 
center {Aa, A5}. The central peak coincident with the stellar position is described by a point 
source with total flux, Fcen, and the offsets of this point source relative to the belt center 
are given by two additional parameters {Aostar, At^star i- The previous ima ging observations 


show that the belt is viewed close to face-on {i = 30°, iGreaves et al.ll2014j) . We £t the SMA 


data directly for an inclination angle, i, and also an orientation on the sky described by a 
position angle, PA, east of north. 


The e Eridani disk spans a large angle on the sky, approaching (or possibly exceeding) 
the half power field of view of the SMA. As a result, the primary beam response has the 
potential to affect the properties derived for the outer regions of the millimeter emission belt. 
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To account for this in the analysis, we multiply each belt model by an accurate, frequency- 
dependent beam model, normalized to unity at the beam center. Appendix A provides a 
detailed discussion of the SMA primary beam shape. For ATCA, with its larger field of view 
at the observed wavelength, the effects of the primary beam shape are much less important, 
and a simple Gaussian provides an adequate description. 


For each set of model parameters, we use the miriad uvmodel task to compute two 
synthetic visibility sets sampled at the same spatial frequencies as our SMA observations, 
corresponding to the two spectrally averaged sidebands (218.9 and 230.9 GHz). The fit qual¬ 
ity is characterized by a likelihood metric, £, determined from the values computed using 
the real and imaginary components at all spatial frequencies (ln£ = —x^/2). This mod- 
eli ng scheme is implernented using the affine-invariant ei isemble MGMG sample r prop osed 
by iGoodman fc Weard (120lOl) and realized in Python by iForeman-Mackev et al.l (120131) . A 
MGMG approach allows us to more effectively characterize the multidimensional parameter 
space of this model and to determine the posterior probability distribution functions for each 
parameter. We assumed uniform priors for all parameters, with reasonable bounds imposed 
to ensure that the model was well-defined: Fbeit > 0 and 0 < < Rout- In addition, we 

constrained the four offset parameters {Aa, A(f, Aostar, ^^star} that describe the belt center 
and stellar position to be within 5" of the offsets predicted from the stellar proper motion; 
this constraint very generously accommodates the uncertainties in the proper motion and 
the absolute astrometry of the observations. 


For the ATGA 7 mm observations, where the e Eridani emission belt is not visible in 
the image, we simplihed the model substantially. We £x all of the model parameters to 
the best-fit values from the analysis of the SMA data, except the total belt flux (Fbeit), the 
total flux for the central peak (Fcen), and two position offsets {Aa, A(5}. In essence, we fix 
the shape of the emission structure to that found from analysis of the SMA data, and we 
determine the 7 mm fluxes for the central component and for the belt component, allowing 
for the possibility of a small positional shift between the SMA and ATGA datasets. We do 
not constrain Fbeit to be positive as for the SMA observations 


3.3. Results of Model Fits 

Table E] lists the resulting best-fit parameter values and their 68% uncertainties deter¬ 
mined from the marginalized posterior probability distributions for both the power law and 
Gaussian belt models fit to the SMA 1.3 mm data. Figure [2] shows a sample of the output 
for the main Gaussian belt parameters, including the marginalized posterior probability dis¬ 
tributions. The la and 2a regions are determined by assuming normally distributed errors. 








- 8 - 


Rc.„ [AU] 



Fig. 2.— A sample of the output from a run of 10^ MCMC trials for the 4 best-fit Gaussian belt 
geometry parameters {Keen, <7, Fheit, and Fcen)- The diagonal panels show the ID histogram for each 
parameter marginalized over all other parameters considered in the model. For each parameter, the 
peak of each histogram is taken to be the best-ht value. The remaining panels show contour plots 
of the 1(7 (red) and 2(7 (gray) regions for each pair of parameters, with the blue crosses marking 
the best-fit values. 


where the probability that a measurement has a distance less than a from the mean value 
is given by erf ■ Both of these functional forms provide good fits to the observed visi¬ 

bilities, with reduced values of about 1.4 (59,924 independent data points, 11 and 10 free 
parameters for the power law and Gaussian models, respectively) for each. Figure |3] shows 
the best-fit models of the 1.3 mm data in the image plane, at the full resolution of the models, 
and imaged like the SMA data, both without noise and with the noise level obtained by the 
observations, which results in patchy outer belt emission very similar to the SMA image in 
Figured] The imaged residuals are also shown in Figure |3l and these are mostly noise (see 
Section [4. 1.51 for further discussion of the residuals from these axisymmetric models). 


A useful way of visualizing the interferometric data and model fits is through the de- 
projected visibi lity function, w hich takes advantage of (near) axisymmetry to reduce the 
dimensionality flLav et al.l 119971) . In particular, the real part of the complex visibilities are 
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Fig. 3.— Images of the best-fit models to the SMA 1.3 mm emission. (Upper): (left) the best-fit 
power law disk model at full resolution, pixel scale ~ Cf!2 (0.8) AU, and (center, left) image of 
the best-fit power law disk model with no noise, (center, right) the best-ht power law disk model 
with simulated noise, and (right) the imaged residuals from the power law disk model, made using 
the same imaging parameters as in Figure [H (Lower): (left) the best-fit Gaussian disk model at 
full resolution, and (center, left) image of the best-fit Gaussian disk model, (center, right) image 
with simulated noise of the best-fit Gaussian disk model, and (right) the imaged residuals from the 
Gaussian disk model, again made using the same imaging parameters as in Figured] The contour 
levels are at 0.4 mJy beam“^ (2 it) intervals in all panels. The ellipse in the lower left corner of the 
residual images indicates the 9'.'2 x T'.'O (FWHM) synthesized beam size. 

averaged in concentric annuli of deprojected {u,v) distance, fZuv, from the center of the 
emission structure. Figure 0] shows this view of the SMA 1.3 mm data together with the 
best-ht power law and Gaussian belt models. The result is a function with a zero-crossing 
null and several subsequent oscillations. Although these SMA observations are missing the 
short (m, v ) spacings needed to sample the peak of the visibility function, the overall shape 
matches nicely that expected for a narrow annulus of emission, plus a small and constant 
positive offset contribution from an unresolved central source. Figure ID also shows the single 
dish hux measurements at the zero-spacing of the deprojected visibility function (with small 
offsets from zero for clarity). It is remarkable that the simple axisymmetric belt models of 
the SMA data appear to provide a very good estimate of the total hux, despite the lack of 
shorter baseline data. This consistency provides indirect support for the basic model of the 
emission distribution. 
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Fig. 4.— The real part of the SMA 1.3 mm visibilities averaged in bins of deprojected (u,v) distance 
from the disk center, compared to the best-fit power law (dashed blue) and Gaussian (solid red) 
disk models plus a centra l point-source. Th e single dish MAMBO-2/1RAM (iLestrade fc Thilliez 


20151 ) and SIMBA fluxes (ISchiitz et al.l 120041 ) extrapolated from 1.2 mm to 1.3 mm are plotted at 


T^uv = 0 kA (slightly offset from zero, for viewing clarity). 


Table m lists the best-fit parameter values from the modeling and their 68% uncertainties 
for the 7 mm data. This model again provides a good fit to the data, with reduced = 1.32 
(151,976 independent data points, 4 free parameters). We note that the best-fit belt flux 
at 7 mm is positive, albeit with low statistical significance; this positive value suggests the 
presence of emission below the detection threshold in individual beams in the image. The 
position offsets between the ATCA and SMA data are small, consistent with zero. 


4. Discussion 

We have performed interferometric observations of the e Eridani system at 1.3 mm and 
7 mm with the SMA and the ATCA, respectively, with baselines that sample to ~ 4" (13 AU) 
resolution. The 1.3 mm image reveals emission from a resolved outer dust belt located ~ 18" 
(60 AU) from the star, and a compact source coincident with the stellar position. The 7 mm 
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image shows only a central peak, detected with greater significance than in the 1.3 mm 
image. We modeled the visibility data assuming two emission components, an outer belt 
with a power law or Gaussian radial surface brightness profile, and a central point source. 
Both functional forms provide good fits to the 1.3 mm data, and we used the best-fit emission 
shape parameters to obtain constraints on the component flux densities from the 7 mm data. 
We now use the new information about these emission components to discuss implications for 
the debris disk and to compare to claims derived from previous millimeter and sub millimeter 
observations with lower angular resolution. 


4.1. Outer Belt 

Several basic properties of the outer emission belt are strongly constrained by the SMA 
1.3 mm observations, including its flux, radial location and width, viewing geometry, and 
departures from axisymmetry. These new constraints bear on the possible presence of unseen 
planets in the system. 


4-1.1. Belt Flux 

The total flux density of the best-fit power law and Gaussian models, is constrained 
to be Fbeit = 16.9^5;g mJy and 17.2~4l4l mJy, respectively. These values are consistent 
with each other, within the uncertainties. They are also consistent with previous single- 



extrapolating the IRAM/MAMBO-2 observation to 1.3 mm gives 14.6 ± 2.9 mJy, while the 
older SIMBA observation gives 18.0 ± 4.3 mJy (see Figure 0]). 


4.1.2. Belt Location and Width 

For the best-fit power law model, the outer radius of the millimeter emission belt is 
determined to be i?out = 80.217;° AU. Previous imaging studies of e Eridani provide esti- 
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mates of the outer radius between 70 — 90 AU, in goo d agreement with this determination 
flGreaves et alJll998l. l2005l. l2014l: iRa.ekman et ahll2009l) . Additionally, we constrain the in¬ 
ner radius of the power law model, = 53. 4l^n AU. This value agrees most closely with 
analysis of 160 pm Herschel observations bv iGreaves et ahl (120141) . which also suggest an 
inner radius of the outer belt of ~ 54 AU. Other single-dish obse rvations indicated that the 
belt extends further inward towards the star, = 35 — 40 AU (IGreaves et al.lllOOSl. 12005 


Backman et al.ll2009l) . The model hts to the SMA 1.3 mm data do not support such a wide 


belt with an inner radius so close to the star. 


The best-£t Gaussian model is characterized by a radial location, i?cen = AU, 

and a width, a = 8.55^3 4 g AU, or FWHM = AR = 2-\/21n(2) x a = 20. 2^^° AU. These pa¬ 
ramet ers are most directly comparable to the belt parameters derived bv iLestrade fc Thilliez 
( 2015 1 using IRAM/MAMBO-2 observations at 1.2 mm with a telescope FWHM beam size 
of lO'.'T; they £t a Gaussian shape to the disk emission radial prohle and obtain a central ra¬ 
dius Rcen = 57±1.3 AU, FWHM = 12" ±1", and infer 8 < AR < 22 AU. This central radius 
is slightly smaller than that derived from the SMA data, but plausibly consistent within the 
mutual uncertainties. Since the lower limit on the width is narrower than implied by the 
£t to the SMA data at the 68% conhdence interval, we examine this potential discrepancy 
more closely. 


The effect of changing the belt location (i?cen) and width (a) is dramatic on the null 
locations in the deprojected visibility function. We expressed our Gaussian model using a 
surface brightness prohle of the form Iu{r) = A x exp[—((r — i?cen)/\/2cr)^]. The Fourier 
transform of a radially symmetric function like this can be expressed by a Hankel transform: 


F{p) = 2%A / R{r)Jo{pr)rdr, 


( 1 ) 


where, p = 271 \/v? -\- = 27iTZuv For a Gaussian ring, there is an exact solution to this inte¬ 

gral involving an inhnite series of hypergeometric functions that can be evaluated numerically 
to hnd the exact locations of the visibility nulls. Fortunately, there is also an app roximate 
solut ion to the Hankel transform using a generalized shift operator (described in iBaddour 
20091 ) that yields the values of the null locations to within 1% of the exact solution. 


F{p) = 271 Aa'^ X exp[-(p\/2(T)V4] x Jo(pi?cen), (2) 

From this simple expression, we can see immediately that for a fixed belt width, decreas¬ 
ing the belt radius (inward towards the star) moves the zero-crossing null locations towards 
larger TZuv For a hxed belt location, decreasing the belt width moves the zero-crossing null 
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locations towards slightly smaller TZ^v (and increases the amplitude of subsequent oscilla¬ 
tions). The best-£t Gaussian model to the SMA data, i?cen = 64 AU, a = 8.6 AU yields 
nulls at TZuv ~ 4, 8, and 14 kA. By comparison, a Gaussian mo del with Rr^-n = 57 AU an d 
Ai? = 8 AU (cr ~ 4 AU), at the lower limit of width suggested by Lestrade fc Thilliez ( 2015 1. 
results in nulls at TZuv ~ 5, 9, and 16 kA, which are signihcantly offset from the data. The 
differences between these ht results may stem from the oversimplihed assumption of a strict 
Gaussian shape for the emission, perhaps exacerbated by deconvolution of the synthesized 
beam resulting from the shift-and-add procedure used to restore the MAMBO map from the 
chopped observations of the e Eridani held. 

The presence of planets can affect the widths of planetesimal belts, through dynamical 
interactions. Given the best-fit Gaussian belt parameters, we can constrain the fractional 
belt width of the e Eridani debris disk to AR/R = 0.3lln'?q. This f ractional width lies 


within the range of 0.1 < AR/R < 0.4 obtained bv iLestrade fc Thilliea (120151) . For com¬ 
parison, the classical Kuiper Belt in our Solar System appears radially conhned between 
40 and 48 AU {AR/R ~ 0.18), hlling the region b etween the 3:2 and 2:1 resonances with 
Neptune, likely the result of its outward migration (lUahn fc Maihotrall2005[l . A narrow ring 
of millimeter emission in the Fomalhaut debris disk (FWHM ~ 16 AU and AR/R ~ 0.1) 
has b een attributed to conhnement by shepherding planets orbiting inside and outside the 
ring (IBolev et al.l 120121) . The best-fit value of the fractional width of the e Eridani belt is 
wider than the Fomalhaut belt and the classical Kuiper Belt, but not as wide as the belt 
surroundi ng the Sun-like st ar HD 107146, which also shows hints of a more complex radial 
structure (IRicci et al.l 120151) . 


The SMA 1.3 mm data do not place any strong constraints on the sharpness of the belt 
edges. This is evidenced by the comparable ht quality for both the sharp-edged power law 
and smoother Gaussian surface density prohles. Sharp edges in the underlying planetesi¬ 
mal distribution might be expected from planetary interactions, as regions within chaotic 
zone boundaries are rapidly cleared. For example, the sharp inner edge of t he Fomalhaut 


debris disk seen in scattered light has long suggested sculpting by a planet (iQuillenI 12006 
Ghiang et al.ll2009l) . In contrast, “self-stirred” debris disks are not expected to show sharp 


edges. In these models, the formation of Pluto-sized bodies initiate collisions that pro pagate 
outward to radii of several tens of AU over Gyr timescales (iKenvon &: Bromlevll2008l ). This 
process tends to produce a radia lly extended planetesimal belt, with an outwardly increasing 
gradient. The specihc models of iKennedv fc WvattI ( 2010l) predict an prohle for the belt 
optical depth. Given the limits of the resolution and sensitivity of the SMA data, model 
htting does not provide a strong constraint on the power law gradient of the belt emission, 
a; = 1.92l°;^t However, if we take this best-ht exponent at face value and assume that 
the emitting dust is in radiative equilibrium with stellar heating, which gives a temperature 
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gradient close to T oc r then this valne implies a rising snrface density prohle, S oc r 
(with large nncertainty). This is similar to the rising snr: 


2.4 


from millimeter observations of the AU Mic debris disk (MacGregor et al.l 120131), as well as 


ace density prohle, S o c fonnd 


the r ising snrface density towards the onter edge of the HD 107146 debris disk flRicci et ah 


20151 ). Since the snrface density of protoplanetary accretion disks decreases with radins, 


this small bnt growing sample of debris disks with rising gradients may point to snpport for 
self-stirred models of collisional excitation. 


4-. 1.3. Belt Viewing Geometry 

In addition to the belt parameters, the data place constraints on the disk viewing 
geometry throngh the inclination and position angle parameters. The position angles from 



ring model to Herschel 160 /rm data and obtain a higher inclination valne, i = 30° ± 5°, 
still compatible with the hts to the SMA 1.3 mm data within the nncertainties. However, 
this difference conld be a sign of backgronnd confnsion affecting the inferences from the far- 
infrared images, or perhaps a wavelength dependence of this parameter dne to the emission 
sampling different grain size popnlations. To assess whether or not snch a higher inclination 
conld affect the other belt parameters in the SMA analysis, we hxed i = 30° and re-ran the 
MCMC model htting procednre; the best-ht parameters are hardly changed (< 3%). 


4 . 1 . 4 . Limits on Belt Stellocentric Offset 


The model hts show no signihcant centroid offset between the belt and central compo¬ 
nent, as might result from the secular perturbations of a planet in an eccentric orbit interior 
to the belt. For planet induced eccentricities, the displacement of the belt centroid from the 
star sh ould be ^ ae, whe re a is the semi-major axis of the belt and e is the forced eccentricity 


(e.g. iChiang et al.ll2009l) . Our modeling allows us to place a robust 3a upper limit on the 
displacement, Arcen 277 = 8.7 AU. Based on a possible far-infrared n orth-south hux asym¬ 
metry attributed to pericenter glow (enhanced emission at periapse), iGreaves et ahl (120141) 
raise the possibility of an additional planet in the e Eridani system with semi-major axis 
within the outer belt a = 16 — 54 AU and eccentricity e ~ 0.03 — 0.3. Given the constraint 
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on the centroid offset from the SMA data, the presence of a giant planet on a wide orbit 
of several lO’s of AU with a large eccentricity, e > 0.1 is disfavored. This is in accord with 


direct imaging constr aints at infrared w’a velengths that preclude planets of about 1 Jupiter 
mass beyond 30 AU fjJanson et al.ll2015l) . However, the effects of a Uranus or Neptune-like 
planet with lower orbital eccentricity is still easily accommodated within the limits. 


4 . 1 . 5 . Limits on Belt Azimuthal Structure 


The azimuthal structure of the e Eridani debris disk has been the subject of much de¬ 
bate. The hrst JCMT/SCUBA 850 /im image of the disk sho wed a non-uniform b rightness 
distribution with several peaks of modest signal-to-noise ratio flGreaves et al.lll998h . Follow¬ 
up JCMT/SCUBA observations from up to 5 years later suggested that three of the peaks 
in the original image appear to move with the stellar prope r motion, and in fac t showed 
tentative evidence of counterclockwise rotation of ~ 1° yr“^ (jCreaves et al.ll2005l) . Several 
other peaks did not appear to mo ve with the star and were presumed to be background fea¬ 


tures. iLestrade fc Thilliezl fl2015l) claim from IRAM/MAMBO-2 observations that the disk 
shows a similar azimuthal structure, with four peaks, in the northeast, southeast, south- 
west, and northwest. H owever, CSO observations at 350 /rm did not conhrm the same peaks 
fiBackman et al.N2009l ). and Herschel observations at 250 /rm (albeit at lower resolution) 
show a relatively smooth emission dist ribution, with the southern end ~ 10% brighter than 
the northern end flGreaves et al.ll2014j) . 


The interest in determining the robustness of the e Eridani clump structure stems from 
the suggestion that the outward migration of a planet could trap planetesimals outside 
of its orbit in mean motion resonances. A variety of numerical simulations show that the 
patt ern of clumps observed in a disk depends on the planet mass and th e resonances involved 


(e.g. iKuchner fc HolmanI l2003l: IWvatt 


2 OO 3 I: IPeller fc Maddisonll2005l) . Thus, the emission 


morphology of a debris belt can be diagnostic of the presence of a planet and its migration 
history. However, other numerical simulations suggest that all azimuthal as ymmetries should 
be ef fectively erased by collisions within debris disks as dense as e Eridani (iKnchner fc Stark 
2010), which would imply that any clumps are spurious, or perhaps background sources. 
Background sources seem particularly problematic to the east of e Eridani in wide held 
Herschel submillimeter images, which likely contributed to confusion at earlier epochs when 
the disk was superimposed on them. 


The SMA 1.3 mm data probes structure at higher angular resolution than the previous 
single dish observations. Moreover, the interferometer naturally provides spatial hltering that 
serves to highlight the presence of any compact emission peaks. After removing azimuthally 
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Fig. 5.— Azimuthal profile of the residual emission after subtracting o ur best-fit Gau s sian m odel 
from the 1.3 mm c o ntinu um emission. The four clumps discussed by I Greaves et al.l (j2005l i and 


Lestrade fc Thilliea (j2015l i are marked by the blue arrows. 


symmetric models from the SMA data, any significant azimuthal structure should be readily 
apparent in the imaged residuals. Figure |5] shows the azimuthal profile of the residual 
image obtained after subtracting the best-fit Gaussian belt model from the data. Each 
point represents the mean brightness calculated in a small annular sector with opening angle 
of 10° and radial range of 10" to 30". Uncertainties are the image rms noise divided by 
the square root of the number of beams in each sector. The only potentially significant 
feature is a peak at 20°, which is visible at the ~ 4a level in the residu al images of Figure [3 


right panels). The loca tions of the previously claimed four clumps (iGreaves et al.l 12005 


Lestrade fc ThilliezI 120151) are marked with arrows in Figure [5l While the azimuthal profile 


of the imaged residuals shows roughly four low signihcance peaks, these are not aligned 
well with the previously claimed clumps, and the separations cannot be readily attributed 
to the previously claimed rotation. However, the si gnal-to-noise of these resid uals is still 
lacking. The signature of the four clumps discussed in lLestrade fc ThilliezI (l2015l) . if they do 
exist, have been weakened in the htting procedure by using an azimuthally uniform ring as 
the model. A more dehnitive statement about low level clumps will require high resolution 
observations with higher sensitivity. 


The imaginary part of the visibilities is very sensitive to the presence of asymmetries in 
the emission structure. Indeed, the effect of the marginally signihcant peak in the residual 
images is apparent. Figure [6] shows the imaginary part of the visibilities binned along the 
M-axis and along the n-axis of the Fourier plane. For a symmetric structure, the imaginary 
part of the visibilities should be zero; however, both of these views show hints of a low 
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Fig. 6.— Complex visibilities of the 1.3 mm emission binned along the t>-axis (top) and the tt-axis 
(bottom). Model visibilities generated by inserting a 1 mJy point source at the position of the 4 cj 
peak seen in the residual map are shown by the blue dashed lines. 


amplitude, sinusoidal oscillation. This structure in the visibilities arises naturally from a 
single, offset emission peak. The Fourier transform of a point source with amplitude. A, 
offset from the center of the image in the u and v directions by xq and yo, respectively, 
is given by G{u,v) = A x exp[—27r(xoM + Vov)] (making use of the shift theorem, whereby 
a shift in the position of a function by a n amount Xn corre sponds to a phase change in 
its Fourier transform by exp(f27rxoM), e.g. Ilsella et al.l (120131) 1: the imaginary part of this 
expression is X{G{u,v)) = Ax sin[27r(a;on + yov)]. That is, a point source offset from (0,0) 
introduces a simple sinusoidal oscillation in the imaginary part of the visibilities. As shown 
in Figure El if we insert a 1 mJy point source at the location of the 4cr peak in the residual 
image [xq ~ 15", y^ ~ 25") into the visibilities, the resulting imaginary visibility curve 
matches very well the shape and scale of the residuals. Thus, this single residual feature 
can account for most of the structure in the imaginary visibilities, or most of the detected 
asymmetry. This feature does not match up with any of the clumps previously identified in 
single dish images. Also, it is positioned clockwise from th e northeast feature, which is not 
consistent with the counterclockwise rotation suggested bv iGreaves et al.l (l2005l) . The exact 
nature of this marginally signihcant peak is uncertain. If it is a background galaxy, then the 
proper motion of the star will move the debris disk away from it, and this should become 
evident in interferometric observations at future epochs. An extragalacti c source should also 


appear compact at the arcsecond level. Recent deep ALMA surveys flHodee et ahl 12013 
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Karim et al.N2013f) have built up statistics describing the number counts of submillimeter 


galaxies expected in a given held of view. At 1.3 mm, the expected numb er of sources wit h 
hux > 1 mJy expected in the 52" primary beam of the SMA is O.dl^g^g f Quo et ahl 12014 ). 
Hence the presence of a background submillimeter galaxy in the image at this hux level is not 
a rare event. Verihcation and characterization of any asymmetric structure in the millimeter 
emission from the belt requires observations with higher resolution and sensitivity. 


4-1.6. Belt Spectral Index 


The long wavelength spectral index of the continuum belt emission rehects the under¬ 
lying dust opacitie s and provides a constraint on the size distribution of the emitting grains 
in the debris disk flRicci et al.l 120121) . which can be related to collisional models. 


For optically thin dust emission, the hux density is given by oc 
where ily(Tdust) is the Planck function at the dust temperature T^ust, oc is the dust 
opacity, expressed as a power law at these long wavelengths, Mdust is the dust mass, and d is 
the distance. At long millimeter wavelengths, for sufficiently high temperatures, the Planck 
function reduces to the Rayleigh-Jeans approximation B,A Tdu<,t) (x w ith api = 2. Thus, 
ttmm ~ «pi + fd, where amm is the millimeter spectral index. iDraind (120061) derived a relation 
between (I, the dust opacity power law index, and q, the grain size distribution parameter: 
fd ^ {q — S)(ds, where (dg = 1.8 ± 0.2 is the dust opacity spectral index in the small particle 
limit for interstellar grain materials, valid for 3 < g < 4 and size distributions that follow 
a power law over a broad enough interval. Combining these relationships provides a simple 
expression for the slope of the grain size distribution, q, as a function of Omm, api, and (dg-. 
q (cTmm Q:pi)/ (dg T 3. 

For e Eridani, the £t to the ATCA 7 mm data places an upper limit on the belt flux 
density of Fbeit < 310 /iJy (3cr). Combining this 7 mm limit with the SMA 1.3 mm measure¬ 
ment provides a long lever arm in wavelength that largely overcomes systematic uncertainties 
and constrains the millimeter spectral index, amm > 2.39. This limit on the spectral index 
results in a limit on the grain size distribution power-law index q > 3.22. 

The derived limit on the grain size distribution power-law index i s consistent wit h 
the classica l pred iction of g = 3.5 for a steady-state collisional cascade (jPohnanvil 119691) . 
Ricci et al.l (120121) obtained a similar result from analysis of the millimeter spectrum of the 
Fomalhaut debris disk, g = 3.48 ± 0.14. This standard collisional cascade assumes that 
collisions in the disk occur between bodies of ident ical tensile strength and velocity dis¬ 
persion, regardless of size. iPan fc Schlichtind (120121) revisited the theory by relaxing the 
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assumption of a single velocity dispersion and solving self-consistently for a size-dependent 
velocity distribution in steady-state. This more complex analysis yields steeper size distri¬ 
butions, with q K, A. Since the e Eridani spectral index is only a lower limit, a steeper grain 
size distribution cannot be ruled out. We note that the best-£t value for the 7 mm flux 
density (Fbeit = HO /uJy) yields q = 3.55 ± 0.36, still shallower than predicted by models 
with size-dependent velocity dispersions. This result, and the more robust measurement for 
Fomalhaut, do not support the steeper size distributions predicted by the collisional models 
with velocity distribution variations. But spectral indices need to be determined for a much 
larger sample of debris disks to draw a dehnitive conclusion. 


4.2. Central Component 


A central source is clearly detected in the SMA 1.3 mm image and ATCA 7 mm image, 
coincident with the position of the star at the time of observation. The source size is below 
the resolution limit of these observations, most clearly evidenced by the lack of fall off at 
long baselines in Figure IS The 1.3 mm flux density of this source is Fcen = 1.081 q4 ® mJy- 
Lestrade fc Thilliez fj2015l ) report a similar value in MAMBO-2/IRAM 1.2 mm data, detect¬ 
ing a central source with flux density 1.2 ±0.3 mJy. These measurements are only marginally 
compatible with expectations for the stellar photosphere at these long wavelengths. The ef¬ 
fective temperature of e Eridani is = 5039 ± 126 K flBaines fc Armstrong! 120121) . and a 
Kurucz stellar atmosphere model (see Backman et ah 2009h predicts a 1.3 mm flux density 
of 0.53 mJy (with 2% uncertainty). The ATCA 7 mm flux density of this same central source 
is Fcen = 66.11 io®5 /iJy, substantially in excess of the stellar photoshere model flux prediction 
of 18 /iJy. The central source persists at a consistent intensity in all of the 8 days of ATCA 
observations, showing no signihcant variability. The mean flux density is 66.5 ± 9.5 /rJy and 
65.9 ± 7.0 /iJy for the June and August observations, respectively. 

In principle, the central 1.3 mm excess emission could be explained by thermal dust 
emission from a warm inner belt. This small 1.3 mm excess, together with Spitzer 24 pm 
and Herschel 70 and 160 pm inner excesses , are consistent with e missi on from a 70 — 100 K 


blackbody, similar to previous inferences bv iBackman et al.l (l2009l) and iGreaves et al.l (120141) 
from the infrared spectrum alone. For reasonable grain sizes, this blackbody emission corre¬ 
sponds to a (very narrow) inner dust belt at 2 — 10 AU, consistent with the size constraint 
from the SMA observations. However, this same blackbody model produces negligible emis¬ 
sion at 7 mm. While previous infrared measurements indicate that there is clearly some warm 
dust present in the system (unresolved in our observations), no inner dust belt scenario can 
also match the substantial 7 mm excess from the ATCA observations. 
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We consider it likely that the unresolved excess emission from the central source arises 
from an additional sfe//ar component, either an ionized corona or chromosphere. The absence 
of variability on the month to month timescale suggests a thermal origin. In particular, the 
millimeter wavelength emission from e Eridani is reminiscent of t he nearby solar-type stars a 
Cen A and B (spectral types G2V and K2V) recently reported bv lhiseau et al.l (1201511 and at - 
tributed to heated plasma, similar to the Sun’s chromosphere. Following Ihiseau et al.l (120131) . 
we calculate the Planck brightness temperature for e Eridani at 1.3 mm and 7 mm, assuming 
the photospheric radius of the star is sufhciently similar at optical and radio wavelengths 
to introduce negligible errors. Optical interf erometry of e Eridani gives a precise measure 
of the stellar radius, i?phot = 0.74 ± O.Oli ?0 flBaines fc Armstrong! 120121) . At 1.3 mm, this 
radius and the excess emission implies Tb = 7800 ± 1400 K, somewhat higher than the op¬ 
tical effective temperature. At 7 mm, however, Tb = 13000 ± 1600 K, very much in excess 
of the photospheric prediction. Indeed, the e Eridani emission follows the sam e trend with 


i ncrea sing wavelength as found for a Cen A and B from ALMA observations. iLiseau et ah 


(120151) measure spectral indices between 0.87 mm and 3.1 mm of 1.62 and 1.61 for a Cen 
A and B, respectively. From the SMA and ATCA data, the spectral index of the central 
component of e Eridani between 1.3 mm and 7 mm is very similar, 1.65±0.23 (where we have 
added in quadrature the ~ 10% flux scale uncertainties at both wavelengths with the la 
errors from model hts). The stellar spectrum clearly starts to deviate strongly from a simple 
optically thick photosphere with a Rayleigh-Jeans spectral index of 2.0, and the contrast 
between the photosphere and the putative chromosphere increases at longer wavelengths. 
While observations of e Eri dani at centi meter wavelengths have s o far provided only upper 
limits, < 80 /rJy at 3.6 cm ( Cudel 1992 1 and < 105 /iJy at 6 cm ( Bower et ah 2009), much 


more sensitive observations are now possible with the upgraded Karl C. Jansky Very Large 
Array, and this would be useful to help constrain the plasma properties. 


5. Conclusions 

We present SMA 1.3 mm and ATCA 7 mm observations of e Eridani, the hrst millimeter 
interferometric observations of this nearby debris disk system, probing to 4" (13 AU) scales. 
These observations resolve the outer dust emission belt surrounding the star, and they reveal 
a compact emission source coincident with the stellar position. We use MCMC techniques 
to ht models of the emission structure directly to the visibility data in order to constrain 
the properties of the two components. The main results are: 

1. The outer belt is located precisely and resolved radially. Caussian and power law emis¬ 
sion prohles each £t the SMA 1.3 mm data comparably well. For the best-£t Caussian 
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model, the belt radial location is Rcen = AU and FWHM = 20.2^g2 AU, 

corresponding to a fractional belt width AR/R = O.Sl^Q^g. This width is at the high 
end of inferences from previous single dish millimeter observations, and wider than the 
classical Kuiper Belt in our Solar System. 

2. The outer belt shows no evidence for signihcant azimuthal structure that might be 
attributed to gravitational sculpting by planets. After subtracting a symmetric model 
from the SMA 1.3 mm data, imaging shows only one low signihcance peak, and its 
location does not correspond to any clumps identihed in previous millimeter and sub- 
millimeter observations of e Eridani. The presence of this feature is consistent with 
source counts for the extraglactic background in the held of view. In addition, the 
SMA 1.3 mm data constrains any centroid offset of the belt from the star to < 9 AU, 
which limits the presence of giant planet perturbers on wide and eccentric orbits in the 
system. 

3. A central source coincident with the star is clearly detected in both the SMA 1.3 mm 
image and the ATCA 7 mm image, and the hux densities of this source exceed ex¬ 
trapolations from shorter wavelengths for the stellar photosphere. While the excess is 
marginal at 1.3 mm, it is highly signihcant- about a factor of three- at 7 mm. The 
stellar spectrum clearly departs from an optically thick photosphere at these long wave¬ 
lengths, with spectral index 1.65 ± 0.23 between 1.3 mm and 7 mm. This spectrum 
cannot be explained by an inner warm dust belt and plausibly results from heated 
plasma in a stellar chromosphere, by analogy with the Sun and a Cen. The high 
brightness temperature at 7 mm of 13000 ± 1600 K for a source of stellar size lends 
additional credence to this conclusion. 

4. Combining the SMA 1.3 mm measurement of the belt hux density with the ATCA 7 mm 
upper limit constrains the spectral index of the emission, amm > 2.39. For conventional 
assumptions about the dust grains, this spectral index corresponds to a limit on the 
slope of the power law grain size distribution in the belt, q > 3.22, consistent with the 
classical prediction of g = 3.5 for a self-similar steady-state collisional cascade. This 
slope is also consistent with the steeper distributions predicted by collisional models 
that allow for size-dependent velocities and strengths. 

These SMA and ATCA millimeter wavelength observations provide the highest resolu¬ 
tion view of the outer dust belt surrounding e Eri at the longest wavelengths to date. But 
deeper observations are still needed to measure radial gradients in the debris disk and to 
reveal substructure due to planets, if present, in order to further constrain scenarios for the 
evolution of planetesimals surrounding this very nearby star. 
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A. Primary Beam Structure of the SMA 

The SMA is composed of eight essentially identical antennas, each 6 meters in diameter. 
The SMA primary beam is thus the power pattern of one antenna. While the primary 
beam shape is often assumed to be a simple Gaussian, the actual shape is determined by 
illumination with a 10 dB taper at the edge of the primary dish, as well as blockage due to 
the secondary mirror. With these considerations, the beam power as a function of offset (in 
arcseconds) from the beam center is given by 


P 


f2Tirx\ /1.840839r\ 
27rr Jo I —^ 1 Jo 


^ Rv 


Rr, 
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(3) 


where, Rp is the radius of the primary dish in meters, Rg is the radius of the secondary dish 
in meters, and x is the offset from the dish center in radians. A complete profile of the beam 
power can be built up using Equation [3] at discrete offset positions. Note that this expression 
does not take into account additional practical factors, such as receiver alignment, pointing 
jitter, and departures from perfect focus, which act to distort the primary beam shape. 

Since the emission extent of the e Eridani disk is comparable to the half power size of 
the SMA primary beam pattern, we constructed a complete beam model for use in our mod¬ 
eling procedure. The FWHM of accurate beam models for the LSB (218.9 GHz) and USB 
(230.9 GHz) are 53"6 and Sff.'S, respectively. For comparison, the FWHM for a uniformly 
illuminated circular aperture antenna is given by 1.22X/Da, where Da is the antenna diam¬ 
eter. For the SMA antennas, this predicts SCf.'d and 47"4, for the LSB and USB, respectively, 
narrower than the FWHM of the accurate beam models. Note that tasks within the miriad 
software package assume a Gaussian beam for the SMA with a FWHM given by a uniformly 
illuminated circular aperture. 
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Table 1. Submillimeter Array Observations of e Eridani 


Observation 

Date 

Array 

Gonfig. 

Projected 
Baselines (m) 

HA 

Range 

225 GHz atm. 
Opacity®" 

2014 July 28 

Subcompact 

6-35 

-3.6,3.3 

0.09 

2014 July 29 

Subcompact 

6-35 

-3.7,2.8 

0.07 

2014 July 30 

Subcompact 

6-35 

-3.5, 7.1 

0.08 

2014 Aug 5 

Subcompact 

6-35 

-3.5,2.2 

0.12 

2014 Nov 19 

Subcompact 

6-56 

-3.6,4.3 

0.07 


Note. — characteristic value for the track measured at the nearby Caltech 
Submillimeter Observatory. The LO frequency for all observations was 225.5 
GHz. 


Table 2. Australia Telescope Compact Array Observations of e Eridani 


Observation 

Date 

Array 

Gonfig. 

Projected 
Baselines (m) 

HA 

Range 

Seeing rms 
(microns)® 

2014 June 25 

H168 

36 - 180 

-3.6,3.6 

150 

2014 June 26 

H168 

36 - 180 

-3.6,3.6 

60 

2014 June 27 

H168 

31 - 180 

-4.0,3.6 

70 

2014 June 28 

H168 

31 - 180 

-4.1,0.0 

250 

2014 Aug 2 

H75 

22 - 84 

-4.7,3.9 

80 

2014 Aug 3 

H75 

22 - 84 

-4.4,4.0 

80 

2014 Aug 4 

H75 

22 - 84 

-4.6,4.0 

150 

2014 Aug 5 

H75 

22 - 84 

-3.8,4.0 

130 


Note. — ^ characteristic value for the track measured by the ATCA 
seeing monitor, an interferometer on a 230 m east-west baseline that 
tracks the 30.48 GHz beacon on the geosta tionary communication s 
satellite, OPTUS-B3, at an elevation of 60°. ( Middelberg et allboOG h 
The LO frequency for all observations was 44 GHz. 
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Table 3. SMA Model Parameters 


Parameter 

Description 

Power Law 

Gaussian 



Best-Fit 

Best-Fit 

Tbelt 

Belt flux density (mJy) 

16.9(-73.9,-5.6) 

17.2(-h5.1,-4.5) 

fr 

^ cen 

Central source flux (mJy) 

1.08(-h0.19,-0.41) 

1.06(-F0.34,-0.34) 

Rin 

Belt inner radius (AU) 

53.4(-F6.1,-5.0) 

- 

Rout 

Belt outer radius (AU) 

80.2(-F3.0,-7.1) 

- 

X 

Belt radial power law index 

1.92(-h0.18,-2.94) 

- 

^cen 

Belt center radius (AU) 

- 

64.4(-h2.4,-3.0) 

a 

Belt width (AU) 

- 

8.55(-H2.54,-3.46) 

i 

Belt inclination (°) 

17.9(-M0.2,-15.3) 

17.3(-hl4.2,-14.2) 

PA 

Belt position angle (°) 

3.42(-h23.3,-23.4) 

1.66(-76.70,-6.70) 

Aa 

R.A. offset of belt center (") 

0.01(-F0.65,-1.03) 

0.02(-F0.80,-0.80) 


Decl. offset of belt center (") 

1.63(-h0.86,-0.86) 

1.63(-h0.70,-1.01) 

AOstar 

R.A. offset of star from belt center (") 

-1.18(-h0.65,-1.40) 

-1.22(-h0.87,-1.26) 

A(5star 

Decl. offset of star from belt center (") 

0.11(-F0.90,-1.29) 

0.11(-hl.l0,-1.10) 


Table 4. ATCA Model Parameters 


Parameter 

Description 

Best-Fit 

Tbelt 

Belt flux density (/rJy) 

110.(-765.,-117.) 

fp 

cen 

Central source flux {fiJy) 

66.1(-h6.9,-10.5) 

Aa 

R.A. offset of belt center (") 

-0.87(-F0.72,-0.87) 

A5 

Decl. offset of belt center (") 

0.38(-F0.86,-0.62) 













